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ABSTRACT 

We present the CXOCY J220132. 8-320144 system, which is composed of an edge-on spiral galaxy 
at z = 0.32 lensing a, z = 3.9 background quasar. Two images of the quasar are seen. The geometry 
of the system is favorable to separate the relative mass contribution of the disk and halo in the inner 
parts of the galaxy. We model the system with one elliptical mass component with the same ellipticity 
as the light distribution and manage to reproduce the quasar image positions and fluxes. We also 
model the system with two mass components, disk and halo. Again, we manage to reproduce the 
quasar image positions and fluxes. However, all models predict at least a third visible image close to 
the disk that is not seen in our images. We speculate that this is most likely due to extinction by the 
disk. We also measure the rotational velocity of the galaxy at 2.7 disk scale radius to be Vc = 130 ± 20 
km s~^ from the [Oil] emission lines. When adding the rotational velocity constraint to the models, 
wc find that the contribution to the rotational velocity of the disk is likely to be equal to or larger 
than the contribution of the halo at this radius. The detection of the third image and a more accurate 
measurement of the rotational velocity would help to set tighter constraints on the mass distribution 
of this edge-on spiral galaxy. 

Subject headings: dark matter — galaxies: spiral, structure, halos - gravitational lensing 



1. INTRODUCTION 

The current picture of galaxies is that they are com- 
posed of baryons (stars, gas, dust) and non baryonic 
dark matter. While observationally the distribution of 
the baryons can be studied, it is difficult to probe how 
the dark matter is distributed compared to the baryons. 
Traditionally, the study of galaxy dynamics has been the 
strongest proof of the existence of dark matter in galaxies 
(e.g., Rubin & Ford 1970; Faber & Gallagher 1979; Rubin 
et al. 1985). But given the inherent degeneracies in the 
inversion of dynamical data to obtain density profiles, it 
is hard to measure how the dark matter is distributed. 

Another approach is to use gravitational lens systems. 
For strongly lensed quasars (QSOs), the geometry and 
photometric properties of the lens system depend on the 
projected mass inside the lensed QSO images and there- 
fore can be used to constrain the mass distribution. 

For the particular case of spiral galaxies, generally 
composed of a bulge, a disk and a dark matter halo, 
the rotational velocity curves cannot disentangle the rel- 
ative contributions of the different mass components in 
the inner (luminous) parts of the galaxy. On the other 
hand, gravitational lenses offer the possibility of doing 
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so, especially in the case of edge-on spirals. 

There are only five spiral gravitational lenses known 
so far (B0218-F357: Patnaik et al. 1993; PKS 1830- 
211: Pramesh Rao & Subrahmanyan 1988; Q2237-I-0305: 
Huchra et al. 1985; B1600-f434: Jackson et al. 1995; 
PMN J2004-1349: Winn et al. 2001). Only one of them 
is seen edge-on, B1600-I-434. However, none of them is 
ideal for studying the relative mass contribution of the 
different components. 

Here, we report the discovery of an edge-on spiral lens 
galaxy as part of the CYDER survey, which may be the 
best spiral lens system to separate the relative mass con- 
tributions of its constituents. The Calan-Yale Deep Ex- 
tragalactic Research (CYDER) survey (Castander et al. 
2003b; Treister et al. 2005, hereafter T05) is an opti- 
cal and near-infrared imaging and spectroscopic program 
carried out in archived, moderately deep Chandra fields. 
CXOCY J220132. 8-320144 (hereafter CY 2201-3201) is 
one of the faint X-ray sources detected by CYDER in 
its Dl field. Optical follow-up of this field seemed to 
suggest that CY 2201-3201 was an edge-on spiral with a 
very bright nucleus. However, optical spectroscopy and 
imaging of this source in good seeing conditions revealed 
its true nature. CY 2201-3201 is a lensing system, show- 
ing two images of a distant z — 3.9 quasar being lensed 
by an edge-on spiral galaxy at z — 0.32. 

The paper is structured as follows. In §2 we present 
in chronological order the observational data gathered 
to characterize this system. In §3 we analyze the ob- 
servational constraints obtained, and discuss them in §4. 
Finally we draw our conclusions in §5. Throughout, we 
assume flo = 0.3, = 0.7 and Ho — TOhjo km s~^ 
Mpc-^ 

2. OBSERVATIONS 
2.1. X-Rays 

CY 2201-32 was observed by the Chandra X-Ray Ob- 
servatory with the ACIS-I instrument for 50.16 ks in 2000 
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Fig. 1. — Top: VLT extracted spectra of image A (upperline) and B (lowerline) of the z=3.90 lensed QSO. For displaying purposes the 
spectrum of component A has been shifted vertically. Bottom: Flux ratio of the two component spectra. In order to avoid contamination 
from the other QSO spectrum, both QSO components extraction apertures are narrow. Notice the small difference in apparent redshift 
due to a slit width considerably wider than the seeing and the two images not being aligned along the slit. This observational configuration 
makes each image of the QSO go through a slightly different part of the grism and therefore the same wavelength is not at the same pixel 
position. 



July as part of the observation of the HCG90 field (PI 
Bothum). We retrieved this image from the archive and 
analyzed it using standard techniques with the CIAO 
package (Castander et al. 2003a, T05). 

CY 2201-32 was detected in the soft band (0.5-2.0 keV) 
with 11.6 ± 3.9 counts and undetected in the hard band 
(2.0-8.0 keV). We computed an X-ray flux of = 1.06 x 
10~^^ erg cm~^ s"""^ and X-ray luminosity of — 1.49 x 
10'''' erg s~i in the soft band (T05). 

2.2. CTIO Optical 

As part of the CYDER survey we imaged the CYDER 
Dl field (T05) with the Cerro Tololo Inter- American Ob- 
servatory (CTIO) 4m telescope using the Mosaic II cam- 
era in 2001 August. We took images for a total of 6600 
and 1800s in the V and / bands, reaching limiting mag- 
nitudes oiV ~ 26.7 and / = 25.1 and effective seeings of 
1.0" and 0.9", respectively. 

CY 2201-3201 appears in these CTIO images as an 
edge-on spiral with a very bright nucleus/bulge with a 
total measured magnitude of Vvega = 21.26 and Ivega = 
19.94.9 

2.3. VLT Spectroscopy 

On 2003 October 31, we took spectra of the X-ray op- 
tical counterparts in the CYDER Dl field with the UT4 
VLT F0RS2 instrument. We used the 300V grism which 

^ Note that these magnitudes differ from the ones reported in 
Treister et al. (2005). The magnitudes presented here are total 
magnitudes while those of Treister et al. (2005) are aperture magni- 
tudes with a relatively small aperture missing considerable amounts 
of flux for this extended source. 



gives a resolution of R ^ 520 (10.5 A) for our 1" slits. 
CY2201-3201 was observed in one of our masks. Given 
the multi-object purpose of the masks, all slits were ori- 
ented north-south. We took five exposures of 1800 s in 
seeing condition of 0.50"-0.75" and bright/grey sky con- 
ditions. We reduced the spectra using standard IRAF^*' 
routines. Surprisingly, for this source we found two spec- 
tra of a QSO at a redshift of z=3.90 in our best seeing 
spectra (see Fig. 1). The two QSO spectra were blended 
in our worse seeing exposures. Our 30 s mask acquisi- 
tion image taken with a seeing of 0.45" confirmed the 
existence of two point sources close to the center of the 
edge-on spiral galaxy, and therefore confirmed the lens- 
ing nature of the system. 

2.4. Magellan Optical Imaging 

We observed the CY 2201-3201 system with the Mag- 
ellan Clay telescope using the MagIC instrument on 2004 
September 8. MagIC has a pixel scale of 0.0691" pixel" 
We took a series of 300 s exposures in three filters (SDSS 
g, r and i) for a total of 1800, 2400, and 1800 s in the g, 
r and i filters respectively. 

We reduced the images using standard procedures in 
IRAF. The resulting combined images have effective see- 
ings of 0.68", 0.65" and 0.62" in 5, r, and i, respectively. 

Figure 2 shows the combined i-band image. The sys- 
tem configuration is composed of an edge-on spiral with 
two images of the QSO at each side of the disk. This is 
one of the expected configurations produced by an edge- 

IRAF is distributed by the National Optical Astronomy Ob- 
servatory which is operated by AURA Inc. under contract with 
the NSF. 
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Fig. 2.— Magellan Magic i-band image of CY 2201-3201. Image 
size is 8" x8". North is up and east is to the left. The QSO images 
are marked A and B, and the galaxy G. 



on spiral galaxy lens (Kccton & Koclianck 1998), except 
that it is missing a third image that should bo located 
close to the disk. 

To obtain tight constraints in the lens modeling it is 
important to know as accurately as possible the position 
and the fluxes of the observed components. In our images 
the seeing is only slightly smaller than the QSO images 
separation. Therefore the overlap between the QSO im- 
ages and the galaxy light distribution makes it hard to 
separate their relative fluxes and determine their posi- 
tions. We have thus used a Monte Carlo Markov Chain 
method (MCMC; Metropohs 1953; Hastings 1970, Gilks 
et al. 1996) to try to determine the positions and fluxes 
(and their errors) of the three components as accurately 
as possible. 

We model the system with three components, two 
point sources (QSO_A and QS0J3) and a spiral galaxy. 
Each point source is modeled with three parameters: the 
a and 6 of its center and its flux. The spiral galaxy is 
modeled as an elliptical exponential with six parameters: 
a and 6 of the galaxy center, flux, axis ratio, position an- 
gle of the major axis and effective radius. Our model thus 
has 12 parameters. 

In order to obtain the best values we generate images 
with varying values of these 12 parameters and convolve 
them with the point spread function as given by bright 
nearby stars. We then compare the resulting model im- 
age with the observed one and compute a goodness of 
fit. We run a MCMC to obtain the best vahies for 11 of 
the 12 parameters. We fix one parameter, the position 
angle of the major axis, because the large elongation of 
the galaxy allows a robust direct determination. Table 1 
presents the parameters used to constrain the lens model 
which result from a combination of the parameters found 
in the three images {g, r, and i). The errors on the fluxes 
have becm artificially increased (see below). 

Table 2 gives the measured magnitudes of the two QSO 
images and the lens galaxy. 



TABLE 1 

Lens Modeling Constraints of CY 2201-3201 

Parameter QSO A QSO B Galaxy 

Ao (arcsec) 0.00 ± 0.02 -0.65 ± 0.03 -0.48 ± 0.03 

A5 (arcsec) 0.00 ± 0.02 -0.51 ± 0.03 0.02 ± 0.03 

Relative flux 1.0 ±0.2 0.9 ± 0.2 

Axis ratio (b/a) ... ... 0.11 ± 0.03 

Position angle (degrees)^ ... ... 140.2 

i?_e (arcsec) ... ... 1.4 ±0.1 

Vc(2.7ro) (km s-i) ... ... 130 ± 20 



Tixed value 

2.5. Magellan Spectroscopy 

On 2004 September 7, 8 and 9, we took spectra of the 

lens galaxy using the Magellan Clay Boiler & Chivens 
spectrograph (B&C). We use a 1" wide slit aligned along 
the major axis of the lens galaxy. We use two grisms: 600 
and 1200 lines mm~^, giving dispersions of 1.6 and 0.8 A 
pixel" ^, respectively. We observed the lens for 30 min- 
utes in 0.8" conditions and 3 hr in 1.2" conditions with 
the lower resolution setting and for 3.5 lirs in 0.7" con- 
ditions with the higher resolution set-up. Unfortunately, 
the slit position moved with respect to the object during 
the observations and the effective on-source exposure is 
shorter. 

We reduced the images using standard procedures in 
IRAF. Figure 3 shows the two dimensional low-resolution 
spectrum at the wavelength of the observed [OH] dou- 
blet. We measure the rotational velocity of the disk using 
the [Oil] doublet. We manage to detect signal out to a 
radius of 3.8" equivalent to 2.7 the effective radius where 
the rotational velocity appears to have already flattened. 
In fact for an exponential mass distribution the maxi- 
mum of the rotational velocity is at approximately twice 
the effective radius (Binney & Tremaine 1987). In or- 
der to obtain a good estimate of the rotational velocity 
we use an MCMC. We model the two-dimensional spec- 
trum assuming an exponential distribution for the galaxy 
light and an arc-tangent function for the rotational ve- 
locity curve. We generate two-dimensional spectra with 
the expression 



I(x, y) = lo exp ( I — Vc arctan ) 

V ro / TT \ Vo J 

where Iq is the normalization, j/o and Xo are the pixel 
coordinates of the center of the galaxy in the two dimen- 
sional spectrum, To is the scale radius of the exponential 
profile, Vo is the "turnover" rotational velocity radius 
and Vc is the asymptotic rotational velocity. We then 
convolve the model image with a Gaussian of the same 
width of the seeing in the spatial direction and with a 
Gaussian of the same width of the spectral resolution in 
the spectral direction. We generate four model images: 
one for the low-resolution data taken on September 7, 
one for the low resolution data taken on September 9, 
one for the high resolution data taken on Septembc;r 7 
and one for the high resolution data taken on Septem- 
ber 8. We compare the model images to the observed 
ones and compute a global of goodness of fit. We run 
an MCMC with two free parameters, Iq and Vc, to ob- 
tain the asymptotic rotational velocity that best fit the 
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TABLE 2 
Magnitudes'' of CY 2201-3201 



Filter 


QSO A 


QSO B 


Galaxy 


9 
r 

ixs 
Ks 


25 06+" " '' 
23.09l°;°t 
22.621°:°^ 
18.75l°;?i 


25 22+"--^^ 
23 20+"-"'' 

22.74i;;:03 

18.96lg;25 


99 „+U.U3 
^^■-^^-0.03 

21.06t°;°i 
17.18«;1 



Tilters g, r, and i: AB system. Ks: Vega system 



observations. Wo fix the rest of the parameters (jjo, Xg 
To and Vq) to the values directly measured from the two 
dimensional spectra (j/o and Xo) or the MagIC images 
(ro). We fix Vq to the same value in pixels of Tq. We ob- 
tain a best value for the asymptotic rotational velocity of 
Vc = 165 ± 25 km s~^. With the arc-tangent rotational 
velocity model used, wc compute a rotational velocity at 
2.7 times the scale radius of 14(2. 7ro) = 130±20 km s"^. 
Most of the signal in the fit comes from data interior to 
this radius and therefore we will use the value at this ra- 
dius to compare to the lens models that best fit the QSO 
image positions and fluxes. 

In addition we measure the galaxy and QSO redshifts 
to be 2 = 0.323 ± 0.001 and z = 3.903 ± 0.002, respec- 
tively. 

2.6. NTT Near-Infrared Imaging 

We observed CY2201-3201 with the ESO New Tech- 
nology Telescope (NTT) on 2004 October 30 in service 
mode. Wc took an exposure of 4185 s in the Ks band. 
The reduced combined image provided by ESO had an 
effective seeing of 0.75" . We calibrated it with standard 
stars. We use the same MCMC method used for the 
Magellan images to calculate the positions and fluxes of 
the three modeled components in this Ks image. The 
calculated positions are much more uncertain than the 
ones calculated with the Magellan images. The fluxes 
nevertheless help us constrain the spectral energy distri- 
butions of the QSO and lens galaxy (see Table 2). 

3. LENS MODELING 

We model the CY 2201-3201 system using the posi- 
tions and fluxes of the two QSO images together with 
the parameters measured for the lens galaxy. We use 
the GRAVLENS software" (Keeton 2001a) developed 
by C. Keeton. The models used are described in Kee- 
ton (2001b). Table 1 summarizes the observational con- 
straints used for the modeling. Note that we have ar- 
tificially increased the error in the measured fluxes to 
allow for micro- lensing and/or QSO variability effects. 
We would like to test the hypothesis that mass is dis- 
tributed similarly to light, and therefore start modeling 
the system with one ellipsoidal mass component forced 
to have the same ellipticity as the one observed in the 
light. We then add a spherical component meant to rep- 
resent the galaxy halo, which we then allow to be slightly 
elongated. We also add external shear to see the role it 
plays in the overall modeling. In this section we discuss 
only the constraints enforced by the observed QSO posi- 
tions and fluxes. We ignore the constraints imposed by 

See http://www.physics.rutgers.edu/~ckeeton/gravlens/ 
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Fig. 3. — Two-dimensional low-resolution B&C Magellan spec- 
trum of CY 2201-3201. The spatial direction runs along the Y-axis 
and the spectral direction along the X-axis. The sky and galaxy 
signals have been subtracted. 

the third and sometimes fourth images predicted but not 
seen and the predicted circular velocity of the model. We 
discuss those constraints in the next section. 

3.1. One Component Elliptical Mass Models 

For simplicity, we start with a singular isothermal el- 
lipsoid (SIE) mass model which gives a flat rotation 
curve. We follow Keeton (2001b) and use his "alpha" 

model. The projected mass surface density is given by: 

k{Q = I C"\ where C = [(1 - + (1 + e)y^Y/^. e 
is related to the axis ratio q hy: = We flx the 

position and ellipticity of the galaxy and only allow the 
normalization to vary to get the best fit to the observed 
QSO positions and fluxes. Figure 4 shows the basic con- 
figuration of the best fit. Three images of the QSO are 
produced at one side of the rotation axis of the galaxy. 
This is a typical spiral galaxy lens configuration (Kee- 
ton & Kochanek 1998), with the two brighter images at 
each side of the disk and the fainter image closer to the 
disk. Figure 4 shows that the position of the QSO in the 
source plane (denoted by the plus sign, +) is very close 
and outside of the radial caustic (solidellipsoidalline) . 
If the QSO were closer to the galaxy center and inside 
this caustic then five images would be produced (with 
the central one highly demagnified). 

Next, we try an elliptical Navarro- Frenk- White (NFW) 
model (Navarro et al. 1996, 1997). The spherical mass 
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Fig. 4. — Fits to the lens configuration (position and fluxes) for one mass component models. Left: SIE model. Middle: NFW model. 
Right: Exponential model. The solid lines show the source plane caustics and the dashed lines the lens plane critical curves. The plus sign 
(+) indicates the position of the QSO in the source plane. The X symbols indicate the two observed positions of the QSO and the dots 
the predicted positions. The numbers at the right and below of the points give the predicted magnification (positive/negative values imply 

positive/negative parity) and the numbers at the left and above the Xs give the re-normalized observed magnifications. 

density in the NFW model is given by p = ^ (i+x)'^ ' 
where x = r/rg. The projected surface mass density 

is K(r) = 2 Kg ^^T-^j whore F{x) is given by equa- 



tion (48) in Keeton (2001b). See also ?; ?; ?; ?; ?. 
Keeton's GRAVLENS code defines an elliptical model 

from this spherical one replacing the polar radius r in 
the projected surface mass density expression by the el- 
lipse coordinate ^ = [x'^+y^/q^Y^^. This model produces 
the same "disk" configuration as the prcvioTis one. We 
fit the QSO positions and fluxes allowing only the nor- 
malization and characteristic radius of the model (rg) to 
vary. We find that a large range of values produce ac- 
ceptable fits. The values of the characteristic radius and 
normalization that produce acceptable fits are degener- 
ate. The trends in this degeneracy are that, when all 
other parameters arc fixed, the larger the characteristic 
radius, the lower the normalization and the smaller the 
radial caustic. That is, the central concentration of the 
model dictates the size of the radial caustic. Figure 4 
shows a typical configuration of this model. 

Finally, we model CY220 1-3201 with an exponential 
mass component (Keeton 2001a, b). The projected sur- 
face mass density of this model is given by /t(^) = 
^ exp{—£^/Rd) where ^ is the same as above, q is the 
axis ratio and Rd is the scale radius of the exponential 
distribution. We fix the scale radius, the position angle 
and the axis ratio to the values obtained with the light 
distribution (Table 1) and solve for the normalization 
that best reproduces the position and fiuxes. The model 
shows the disk configuration as before and reproduces 
the position and fluxes of the two QSO images fairly well 
(see Fig. 4). 

Given that the three models produce acceptable fits, 
one would like to find another constraint to differentiate 
between them. Comparing the three models, we find that 
the position and relative magnification of the third image 
depend on the concentration of the mass models. For the 
most concentrated model (SIE), the position of the third 
image is further away from the galaxy center and the 
relative magnification compared to the brightest QSO 
image is the lowest (30%). For the least concentrated 
model (exponential), the position is closer to the galaxy 



center and the relative magnification is the largest (59%). 

3.2. Two Components Mass Models (Adding a Spherical 

Halo) 

We start with a simple two component model: one 

SIE, as in the previous section, to which we add a sin- 
gle isothermal sphere (SIS). The SIS model uses the al- 
pha model of Keeton (2001b) (the same as the SIE, 
see above) in which e = 0, the axis ratio q = 1, and 
C = [x^ + y^]^/^. We vary the normalizations of the SIS 
and the SIE models, fixing all other parameters, to find 
the best fitting model for the observed data. A wide 
range of normalizations in the models produce accept- 
able fits to the position and fluxes. All acceptable fits 
produce five images of the QSO (the central one strongly 
de-magnified). Only models in which the rotational ve- 
locity of the SIS is approximately 50% larger than the 
rotational velocity of the SIE do not give acceptable fits. 
Figure 5 shows an example of an acceptable fit in which 
both components have similar rotational velocities. Best 
fits are produced for models with larger SIE contribu- 
tion. In general when the contribution from the spherical 
component (SIS) gets larger in the fits, then the follow- 
ing trends are observed: the area enclosed by the radial 
caustic increases and the area enclosed by the tangential 
caustic decreases; the magnification of the QSO images 
gets larger and the relative magnification of the fourth 
brightest image compared to the brightest gets larger. 
On the other hand, the relative magnification of the third 
image compared to the brightest is almost independent 
(staying around 30%) of the ratio of the SIE and SIS 
contributions. The fourth and fifth images disappear at 
almost the point where the SIS contribution is negligible 
and the model is dominated by the SIE component. 

We next model the system with an expectedly more 
realistic mass distribution with two components: an ex- 
ponential distribution (as in §3.1) to account for the disk 
and a spherical NFW (same as above but with the axis 
ratio q = I) representing the halo. We search for models 
that best fit the data, varying the two normalizations and 
the characteristic scale rg of the NFW model, and fixing 
all other parameters. Again, acceptable fits to the posi- 
tions and fluxes are produced for a wide range of model 
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Fig. 5. — Fits to the lens configuration (position and fluxes) for two-mass component models. Left: SIE+SIS. Right: Exponential+NFW. 
The symbols are the same as in Fig. 4. 



parameters. The only models that are rejected are those 
in which the spherical NEW component dominates. This 
occurs when the rotational velocity of the NFW is equal 
to or 25% larger (depending on the concentration pa- 
rameter of the NFW) than the rotational velocity of the 
exponential component at 2.7 disk scale radii. These re- 
jected models produce five images, with the fifth central 
image strongly dcmagnified. The only exception is when 
the combination of the two mass components is such that 
the QSO needs to be close to the caustic to fit the ob- 
served position and fiuxes. In this case the fifth image 
moves away from the center towards the fourth image in 
the image plane. The trends with respect to the rela- 
tive strengths of the exponential and NFW components 
are the following: The larger the NFW spherical com- 
ponent, the worse the fit to the positions and the larger 
the magnification of the QSO images. The relative mag- 
nification of the third brightest image compared to the 
brightest image goes from ^50% when the NFW domi- 
nates to ^59% when the exponential dominates. 

If we model the halos with some ellipticity, the results 
present general trends that are between the previous one 
mass component and two mass components cases dis- 
cussed. The most significant difference in the case of an 
exponential plus elliptical NFW model is the magnifica- 
tion of the resulting images. The total magnification of 
the QSO images is in general lower and the relative mag- 
nification of the third image compared to the brightest 
becomes larger than in the previous cases. 

3.3. Adding the Rotational Velocity Constraint 

In the previous section we restricted ourselves to mod- 
eling the system taking into account the positions and 
fiuxes of the two observed QSO images as constraints. 
We have also measured the rotational velocity of the lens 
galaxy and can use it as an additional constraint. 



For each model we compute the rotational velocity pre- 
dicted at 2.7 To, where To is the scale radius of the ex- 
ponential distribution of the galaxy light. This is the 
value of the radius for which we can measure the rota- 
tional velocity with our spectroscopic data (see §2.5). All 
rotational velocity comparisons are performed at this ra- 
dius, neglecting the rotational velocity dependence with 
radius. 

For the one mass component models, the best fitting 
SIE model to the QSO image positions and fluxes pre- 
dicts a lower rotational velocity than observed, but con- 
sistent within the 1 a error. The exponential and ellipti- 
cal NFW models predict larger rotational velocities than 
observed, which are consistent only at the 2 a level. 

If we include the rotational velocity in the overall min- 
imization of the two component models, we find that for 
the SIE-I-SIS model the best fit occurs when the rota- 
tional velocities of the SIS is approximately 20% lower 
than the rotational velocity of the SIE. Unacceptable 
models happen when the rotational velocity of the SIS is 
~ 5% larger than the rotational velocity of the SIE. For 
the exponential-|-NFW case, we do not find acceptable 
models that satisfy the three constraints: positions and 
fluxes of the QSO images and rotational velocity of the 
galaxy. The general trend is that the larger the relative 
contribution of the spherical component (halo), the lower 
the rotational velocity predicted to fit the positions and 
fluxes. 

4. DISCUSSION 

CY 2201-3201 is an edge-on galaxy lens. Such systems 
offer the possibility of decomposing the relative mass con- 
tributions of the disk, bulge and halo. In our optical im- 
ages CY 2201-3201 appears as a bulgeless edge-on spiral 
galaxy, producing two images of a background z = 3.9 
QSO. We have therefore modeled the system with one 
(disk) and two (disk-l-halo) mass components, neglecting 
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the bulge. 

All viable models explored produce cither three or four 
visible images. However, our optical images only show 
two images of the QSO (Fig 2). We observe the two 
images that are farther away from the disk but miss 
the third (and in some configurations the fourth) im- 
age located very close to the disk. We have investi- 
gated whether dust extinction could be responsible for 
this missing image. The colors of the QSO images A 
and B (Table 2) are in fact redder than expected for a 
typical QSO at that redshift. Although low in signal- 
to-noise ratioand covering a small wavelength range, the 
QSO spectrum does not reveal indication of strong in- 
trinsic absorption in the QSO itself. One is then led 
to conclude that the redder colors are due to dust ex- 
tinction between the QSO and the observer. The most 
probable source of attenuation comes from the disk itself. 
We have computed the amount of dust extinction at the 
galaxy redshift {z = 0.323) necessary to explain the ob- 
served colors assuming a mean QSO spectrum at this 
redshift, no intrinsic absorption at the QSO itself and 
the value of the local Galactic extinction Ay = 0.092 
(Schlegel et al 1998). We find that the QSO image A 
requires an extinction by the galaxy disk of Ay = 1.65 
and the image B an extinction of Ay = 1.55. Images 
A and B are seen at a projected distance of 0.35" and 
0.47" (1.7 and 2.2 hj^^ kpc) of the disk plane respec- 
tively. The inferred values of the extinction by the disk 
are then likely taking into account thee projected dis- 
tances to the disk. Our i band image is the one with the 
highest signal-to-noise detection of the QSO images. We 
have inserted a third image of the QSO in our model im- 
ages at the position given by our lens models^'' and verify 
through our MCMC method the maximum flux it could 
have without being significantly detected. If the ratio of 
expected flux to maximum observed flux are completely 
due to absorption by dust at the lens galaxy, then there 
are at least 3 mag of extinction more at the position of 
the third image than at the brightest. Taking into ac- 
count the possible extinction at the brightest image, the 
third image would suff'er Ay > 4.7 of extinction going 
through the lens galaxy. 

In our models of the lens, we have first explored the 
assumption that the mass distribution has the same el- 
lipticity as the light. We have tried three one- mass com- 
ponent models fixing the galaxy center, axis ratio and 
position angle to those observed in the optical. We have 
fit an SIE, an exponential and an NFW model. All 
three models reproduce the two QSO images and predict 
a third image that is unobserved. The predicted rota- 
tional velocities are consistent with the observed value 
at the 1 a level for the SIE model and at the 2 cr level for 
the exponential and elliptical NFW models . We have 
also tried other methods to estimate the galaxy mass. 
We have fit the observed galaxy photometry with the 
PEGASE synthesis evolutionary models (Fioc & Rocca- 
Volmerange 1997). We obtain an absolute magnitude of 

^'^ In fact, there are some possible configurations where the fifth 
central image, normally largely de- magnified, should also be visible. 

Note that the exact position depends on the model. In fax;t it 
can be used to distinguish between them. Nevertheless for the pur- 
pose used here, the position can be considered similar throughout 
all plausible models. 



Mr = —20.37 + 5 logh7Q, which is a factor 3 fainter 
than L* (Blanton et al. 2003). The stellar mass-to- light 
ratio(M/L) of the best fitting model is M/Ly ~ 4.0, so 
its stellar mass would be M ~ 3.8 x 10^° /if^f Mq. The 
expected rotational velocity for such a mass at 2.7 times 
the scale radius is Vc = 60 — 85 km s^^ (depending on the 
mass model). The stellar mass by itself is thus insufficient 
to produce the observed (or predicted) rotational veloc- 
ity. If we assume the local Tully-Fisher relation (Tully & 
Fisher 1977; Pierce & Tully 1992; Bell & de Jong 2001; 
Verheijen 2001) neglecting evolution and the measured 
absolute magnitude, we obtain a value for the expected 
rotational velocity of Vc ~ 125 — 150 km s~^, which is 
consistent with the measured rotational velocity. 

We have also studied more realistic models. The galaxy 
does not appear to have a bulge in the optical images, 
and therefore we have modelled the system with two mass 
components: one for the disk and one for the halo. We 
have tried an SIE-I-SIS model and an exponential-|-NFW 
model. These models reproduce the positions and fluxes 
of the two observed QSO images. However, the SIE+SIS 
model predicts four visible images and the best fitting 
exponential-|-NFW predicts three visible images. None of 
these predicted additional images are seen in our images 
(see above). 

We measure a rotational velocity of Vc = 130 ± 20 km 
s""'^ at 2.7 disk scale radius for the lens galaxy. The 
SIE-I-SIS model predicts values of the rotational veloc- 
ity consistent with this value for certain combinations of 
the relative contributions of the SIE and SIS components 
(see §3.3). However, the exponential-!- NFW model pre- 
dicts higher values of the rotational velocity if we fit the 
observed position and fluxes of the QSO images. Our 
Magellan spectroscopic data were obtained on three dif- 
ferent nights with two different settings. As explained in 
§2.5, we have four sets of same night/grating data. The 
individual fits to the rotational velocity at 2.7ro for each 
set of data are 140 ± 40, 160 ± 45,90 ± 40, and 130 ± 35 
(low-resolution September 7, low-resolution September 9, 
high-resolution September 7 and high-resolution Septem- 
ber 8, respectively). The relative dispersion of these val- 
ues may hint at a possible underestimation of the errors. 
It this were the case, the range of allowed Vc values would 
be larger and the exponential-|-NFW model prediction 
would still be viable. 

We have also explored other eff'ects that can aff'ect our 
modeling. Our system is likely to be influenced by some 
external shear which will contribute to the image separa- 
tion but not to the rotational velocity of the lens galaxy. 
In fact, CY2201-3201 lies 7' away from the Hickson com- 
pact group HCG 90. We have computed what would be 
the external shear produced by the group assuming it is 
modeled with a SIE with the same velocity dispersion 
as measured from the galaxy members. HCG 90 is rel- 
atively small and the external shear induced in the CY 
2201-3201 system is negligible for our purposes. Apart 
from HCG 90, CY 2201-3201 appears to be isolated and 
not in any group, cluster, or obvious large-scale struc- 
ture. In fact, as part of the CYDER survey we have 
obtained spectra of several sources in the field and, with 
the limited spectroscopic data we have, not found any 

Taking into account only the K-correction. Mr = —20.11 -|- 
5 loghjo adding an evolutionary correction. 
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sign of a massive structure. 

Another possible flaw in our modeling could be that 
the galaxy center is miscalculated. If the galaxy cen- 
ter is much closer to the QSO images than the position 
we have measured then most of the discussed configura- 
tions would no longer apply and the system would dis- 
play other image configurations, which can in fact place 
tighter constraints on the relative contribution of the halo 
and disk to the total mass budget in the central regions 
of the lens galaxy. However, the seeing and pixel size of 
our images and the consistency of the galaxy center in 
our different filter images make us believe that the true 
center if different from the one measured should not be 
very far off. 

5. CONCLUSIONS 

Wc have presented the discovery and subsequent follow 
up observations of the CY220 1-3201 system composed of 
an edge-on spiral at ^; = 0.323 splitting a background 
z = 3.903 QSO into two observed images each at opposite 
sides of the disk. 

Wc have modeled the system with one (disk) and two 
(disk+halo) mass components. The most likely config- 
uration is the "disk" configuration with three images of 
the QSO one at each side of the disk and the fainter one 
approximately behind the disk. There are also possible 
configurations that produce four or even five (this one 
unlikely) observable images. However we only observe 
two images. We have discussed the possibility that the 
third (and fourth, if existent) image is extincted by the 
disk. We estimate that an Ay > 4.7 at th epredicted 
position of the third image at the galaxy lens redshift is 
required to be consistent with our observations. 

We have measured the rotational velocity of the lens 
galax;y to be 14 = 130 ± 20 km s^^ at a radius 2.7 times 
the scale radius of the galaxy exponential light distribu- 
tion. If wc use an SIE+SIS model to fit the QSO image 
positions and fluxes and this value of the rotational veloc- 
ity, we find that the contribution by the SIE (disk) to the 
Vc at this radius is required to be the same or larger than 
the contribution of the SIS (halo). If we use an exponen- 



tial+NFW model then we are unable to reproduce this 
value of the rotational velocity if we fit the positions and 
fluxes of the QSO images. We have speculated whether 
we have underestimated the error in our determination of 
the Vc which would make the exponential-|-NFW model 
viable. 

CY 2201-3201 is the best Icnsing spiral galaxy known 
to date that can be used to disentangle the contributions 
of its different mass components. Unfortunately our cur- 
rent follow up observations are not constraining enough 
to elucidate between different possible models. More ac- 
curate image source positions, the discovery or not of the 
predicted third image and the measurement of its prop- 
erties and a precise rotational velocity would make this 
system fulfill its potential. 

CY 2201-3201 has been awarded Hubble Space Tele- 
scope (HST) ACS time in cycle 13. Wc have also been 
awarded more spectroscopic time at Magellan. We ex- 
pect that the higher quality images and extra spectra 
will help us improve our modeling and place strong con- 
straints on the relative contribution of the disk and halo 
mass components. 
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